Edge states in Graphene: from gapped flat band to gapless chiral modes 
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We study edge-states in graphene systems where a bulk energy gap is opened by inversion symmetry break- 
ing. We find that the edge-bands dispersion can be controlled by potentials applied on the boundary with unit 
cell length scale. Under certain boundary potentials, gapless edge-states with valley-dependent velocity are 
found, exactly analogous to the spin-dependent gapless chiral edge-states in quantum spin Hall systems. The 
connection of the edge-states to bulk topological properties is revealed. 
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One of the most intriguing phenomena in solid state physics 
is the existence of edge-state on the boundary of a 2D system. 
The edge-state can have distinct properties from the bulk band 
and play important roles in transport. When the system has a 
band structure in which empty bands and fully occupied bands 
are separated by an energy gap, cuiTent can not flow in the 
bulk. However, this does not dictate the system to be a simple 
insulator, as conduction may still be allowed by edge-states 
on the boundary. Well known examples are the quantum Hall 
effect (QHE) and the quantum spin Hall effect (QSHE), where 
gapless chiral edge-states are robust channels with quantized 
conductance fH [j, S B H Si ■ On the other hand, the property 
of the edge-states are intimately related to the property of the 
bulk band. For example, in QHE and QSHE, the existence 
and chiral nature of gapless edge-states are found to be the 
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consequence of bulk topological orders 17. 8ll9l ll0lll 1LI12I1 . 

The recent realization of graphene in laboratories has at- 
tracted extensive interests to this 2D lattice HEIQ. Its 
bulk band structure has two degenerate and inequivalent val- 
leys centered at the corners of the first Brillouin zone (known 
as the Dirac points). The free standing graphene crystallite is 
a zerogap semiconductor where the conduction and valance 
band touch each other at the Dirac points. A bulk energy 
gap can be opened by breaking the inversion symmetry, e.g. 
by a staggered sublattice potential in single layer graphene, 
or by an interlayer bias in graphene bilayer The neutral 
graphene system then becomes a normal insulator, with van- 
ishing Hall conductance and spin Hall conductance. Nev- 
ertheless, it does acquire a nontrivial bulk topological prop- 
erty: the two inequivalent valleys each caiTy a non-zero topo- 
logical charge with opposite signs, giving rise to the valley- 
dependent Hall effect il6l IitII . Not by coincidence, edge- 
states also exist in such graphene systems with peculiar behav- 
iors |18, 19, 20, 21, 22, 23, 24]. For graphene sheet terminated 
with zigzag edges, edge-states form a one dimensional band 
which connects the two Dirac points with completely flat dis- 
persion. Because of its unusual dispersion relation, the edge- 
bands in graphene have been exploited for a number of in- 
teresting phenomena including valley-filtered transport ll25ll . 
magnetism |26, 27, 28, 29], and superconductivity 13011 . 

In this work, we study edge-state behaviors in graphene 
systems with broken bulk inversion symmetry. We find that 




FIG. 1: (Color online) (a) Schematic illustration of a single layer 
graphene with zigzag edge on the left boundary. The two triangular 
sublattices are denoted by • and o respectively. A unit cell contains 
one atom from each sublattice (denoted by the oval), (b) Band struc- 
ture of graphene with a zigzag edge, (c) Band structure of graphene 
with a bearded edge (i.e. with column L\ removed in (a)). In (b) and 
(c), only the edge-states on the left boundary are shown (red curves). 
A staggered sublattice potential A/2 = 0.1 is universally applied in 
the bulk and on the boundary (see text). 



the edge-band dispersion can be continuously changed by 
simply tuning the on-site energies on the boundary of the sys- 
tem. Under certain values of the boundary potential, the edge- 
band can either completely merge into the bulk band structure, 
or become valley-dependent gapless chiral modes. In the lat- 
ter case, for single layer graphene (bilayer graphene), each 
boundary of the system carries two (four) gapless edge-states, 
one (two) at each valley, with opposite velocities. Therefore, 
when the Fermi energy lies in the bulk gap, we have quan- 
tized valley cuiTent flow on the boundary provided that inter- 
valley scattering can be neglected. This is in close analogy 
to the quantized edge spin current in QSHE jH 0, S. Inter- 
valley scattering in graphene is well suppressed by the large 
momentum separation 11311 l32l 13311 . While these graphene 
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systems do not have nontrivial Z2 topological order in the 
bulk g M, 

[m [Till , the edge-states can nevertheless have 
the same features as those edge-states in QSHE, with the val- 
ley index playing the role of the spin. We further prove that 
these gapless edge-states are equivalent to the gapless chi- 
ral modes confined by topological domain wall in graphene, 
whose existence and chiral nature have been related to the 
valley-dependent bulk topological charge ll34l 351. There- 
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fore, the peculiar edge-states on the boundary indeed share 
the same topological origin as the valley-dependent Hall ef- 
fect in the bulk. This also opens up a new perspective to the 
existence of gapless chiral edge-states by contrasted topolog- 
ical charges in different regions of the Brillouin zone. 

In the tight binding approximation with nearest neighbor 
hopping energy t, a single layer of graphene can be described 
by the following Hamiltonian 



(1) 



where ^Yli{i j) sums over only nearest neighbor pairs. We con- 
sider a general situation where the bulk lattice can be subject 
to a staggered sublattice potential: Ui = A/2 for sublattice 
and Ui = —A/2 for sublattice o (see Fig. [Tfa)). On the 
boundary, graphene sheet terminated with the zigzag or the 
bearded edge preserves the two-valley band structure, and the 
edge-states form well known flat ID bands which connect the 
two Dirac points [18, 19, 20]. We find that the edge-states 
are also gapped when the staggered sublattice potential is uni- 
versally applied to all sites including those on the boundary. 
In Fig. [Tfb) and (c), for clarity, only the edge-states on the 
left boundary are shown, while we note that an edge-state 
on the left boundary with wavevector ky and energy E al- 
ways has a counter part on the right boundary with ky and 
—E. For the zigzag edge, the flat band appears in the region 
ky G [|7r, |7r], while for bearded edge, it appears in the com- 
plementary region ky ^ [|7r, Itt]. Here and hereafter, we nor- 
malize all length scale by the lattice constant a and all energy 
scale by the hopping energy t. Below, we focus on systems 
with zigzag edges, while qualitatively the same behaviors are 
always found in the systems with bearded edges. 



Although the edge-states in the region ky G 



[fTT, Itt] all 

have the same energy (Fig. [TJb)), their degree of localiza- 
tion in X direction varies. The edge-states near ky ^ n are 
almost completely localized on the outermost carbon atoms 
while those near the two Dirac points are much more spread 
into the bulk. Thereby, the energy response of the edge-states 
to potentials applied with unit cell length scale on the bound- 
ary will be different, which forms the basis for controlling the 
dispersion of the edge-band. The edge-bands on the two op- 
posite boundaries can be individually controlled, and in Fig.|2] 
we demonstrate this controllability by tuning the on-site en- 
ergy U{Li) on the outermost column Li (cf. Fig. [Ha)). Since 
the edge-state at ky = tt is almost completely localized on 
column Li, its energy is just given by U{Li). Therefore, the 
edge-band dispersion bends upward when U{Li) increases 
[Fig. |2a-d)]. When [/(Li) — 1, the edge-band completely 
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FIG. 2: (Color online) Band structure of single layer graphene with 
zigzag edges. A staggered sublattice potential A/2 = 0.1 is applied 
throughout the sheet, while the on-site energies U{Li) of the out- 
ermost column Li (see Fig.[TJa)) are tuned to different values from 
positive to negative, (e) with U{Li) — A/2 — 0.1 is the same plot 
as Fig. [lib). The inset in (g) is a blow up of the bulk gap region near 
ky — iyr. Only the edge-states on the left boundary are shown (red 
curves). 



merges into the bulk conduction band continuum. With fur- 
ther increase of U{Li), the edge-band reappears on top of 
the bulk conduction continuum. Meantime, a new edge-band 
starts to peal off from the bulk valance continuum in the com- 
plementary region ky ^ [|7r, |7r] (Fig.|2];b)). This new edge- 
band will approach a flat dispersion at sufficiently large posi- 
tive U{Li) (Fig.|2ja)), which is expected since carbon atoms 
in Li column is then effectively decoupled from neighboring 
columns and the graphene sheet effectively terminates with 
bearded edge (cf. Fig. [lie)). The situation is similar when 
U{Li) is decreased to negative values [Fig.|2|f-i)]. The edge- 
band dispersion bends down initially and traverses the gap 
when[/(ii) < —0.1. AtU{Li) = —1, the edge-band merges 
into the bulk valance continuum. Most significantly, two gap- 
less edge-modes with opposite velocity appear in the vicinity 
of the two Dirac points respectively (see Fig.|2g) and inset). 



The gapless edge-states usually reflect non-trivial topologi- 
cal orders in the bulk. In QHE, it is well known that the num- 
ber of gapless chiral edge-states is given by the bulk Chern 
invariant [8]. This connection has been generalized to QSHE 
recently [l^ IJI, where the existence of gapless chiral 
edge-states are dictated by a generalized definition of bulk 
Chern invariant [9], also known as the Z2 topological invari- 
ant ifioi [TTh . To understand the edge-states in the present 
graphene system, we shall look into its bulk properties. 

In the bulk, we can make a Fourier transform Cj ,/o = 
J2k'^k,»/o^^'^^' with R4 being the lattice vector, and the 
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FIG. 3: (Color online) Schematic illustration of a topological domain 
wall. • indicates carbon atoms with higher on-site energy -L^, and o 
indicates carbon atoms with lower on-site energy — . The index of 
lattice sites k, . . . are arranged symmetrically about the domain 
wall in the regions to the left and to the right (denoted as jC and TZ 
respectively). 



Hamiltonian in Eq. ([T]i becomes 



A/2 V{k) 
V*ik) -A/2 



(2) 



where V{k) = - (l + e"*'" "! + e^^'^ ^^), and ai,2 are the 
two primitive translation vectors (see Fig. [11^)). Eq. (|2|i has 
the solutions of a conduction band \uc,k) and a valance band 
|u„ fc). In order to describe the bulk topological property for 
the graphene system as an insulator, it is convenient to intro- 
duce a gauge potential and a gauge field in the valance band, 
defined as A(fc) = (M„.fc|iVfc|u^.fc) and J7(fc) = VfcX ^re- 
spectively. This gauge field, known as the Berry curvature, is 
analogous to a magnetic field in the crystal momentum space. 
Its integral over a /c-space area yields the Berry phase of an 
electron adiabatically going around the boundary of the area, 
which is similar to the relationship between a magnetic field 
and the Arharonov-Bohm phase. In 2D system, the Berry cur- 
vature vector is pointing out-of-plane and the Chern invariant 
is given by the flux of the Berry curvature threading the en- 
tire Brillouin zone C — Jg^ (Pkn{k). For the graphene 
system described by Eq. (|2]l, one finds the Berry curvature 
has a distribution sharply centered at the two Dirac points, 
ri(q) = Tz 2{A^+tq2)3/2 ^ where = ± is the index for the 
two valleys and q is the wavevector measured from the Dirac 
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points 11161 llTT . We find the two valleys each carry a topo- 
logical charge oi N3 = ^ J (Pq^l{q) — ir2Sgn(A), corre- 
sponding to a half-quantized valley Hall conductivity in the 
bulk m 



1711 . However, the Chern invariant is zero because 
of the time reversal symmetry in the system. It can be further 
shown that the Z2 topological invariant also vanishes when the 
bulk gap is from inversion symmetry breaking only. Clearly, 
the gapless edge-states we found here do not have the same 
origin as those in the QHE and QSHE. 

In searching for the topological origin of these edge- 
states, we notice that similar gapless valley-dependent chiral 
zero-modes also exist at topological domain walls in biased 
graphene bilayer, and the number and chiral nature of these 



modes have been related to the valley-dependent topological 
charge Nj, in the bulk 1 34] . More specifically, in each valley, 
the topological charge v of the zero-modes, i.e. the number 
of zero-modes moving in +j/ direction minus the number of 
zero-modes moving in ~y direction, is equal to the difference 
of the bulk topological charges of the vacua on the two 
sides of the interface, v = iV3(right) - N3(lcft) |34,'35]. For 
the single layer graphene under staggered sublattice potential, 
we can similarly consider a topological domain wall repre- 
sented by a sharp kink in the order parameter A as shown in 
Fig. [3] Gapless chiral modes with topological charge v — 
is thus expected in the domain wall region ||36|] . Interestingly, 
V ~ Tz 'v& exactly the topological charge of the gapless edge- 
states shown in Fig.[2tg), when the on-site energy of the atoms 
on the outermost column is equal to the nearest neighbor hop- 
ping energy. Below, we show the intrinsic connection between 
the gapless zero-modes in the domain wall and the gapless 
edge-states on the boundary. 

Graphene with the above mentioned topological domain 
wall can be described by the following Hamiltonian 

a = - XI + ^ C/,,sc[sCi,s (3) 

where the index i,j here run over lattice sites on one side 
of the domain wall. Cj^A = (cj.Tj — Cj^c) and Cj.s = 
(cj.K + Cj,£ ) describe respectively the antisymmetric and 
symmetric combination of wavefunction on the two sites 
(i,£) and («,7?.) located symmetrically on the two sides of 
the domain wall (see Fig. O. For i in the nearest column to 
the domain wall (Li in Fig.O, we have Ui^s = — 1 + and 
Ui_A = 1 + respectively, while for all other lattice sites, 
both Ui^A and J7i 5 are equal to on sublattice • and equal 

to — on sublattice o. In this way, we establish the equiva- 
lence between an extended graphene with a domain wall and 
a semi-infinite graphene with one zigzag edge subjected to 
certain boundary potential. In particular, the gapless chiral 
zero-modes in the domain wall are of symmetric wavefunction 
and are equivalent to the gapless edge-states on the boundary 
when on-site energy of the outermost column is — 1 + (cf. 
Fig. |2|g)). Therefore, the number of gapless edge-states in 
the latter situation is also determined by the bulk topological 
charge 7V3 in the two valleys. 

We note that the edge-states can evolve continuously from 
the gapless chiral modes near the Dirac points to the flat 
dispersion band connecting the two Dirac points by tun- 
ing the magnitude of the boundary potential only, during 
which the bulk property is unchanged. Therefore, these edge- 
states spectra shall all have the same origin from the valley- 
dependent bulk topological charge, and in particular, the num- 
ber of edge-band is determined by the value of N-^. 

Finally, we turn our attention to a different graphene sys- 
tem, the bilayer graphene with Bernal stacking. By applying 
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termined by the valley-dependent topological charges in the 
bulk. 

The work was supported by NSF, DOE, the Welch Founda- 
tion, and NSFC. 



FIG. 4: (Color online) (a) Schematic illustration of a graphene bi- 
layer. (b) Band structure of a biased graphene bilayer with zigzag 
edges. The on-site energy is universally —0.015 in the upper 
layer and 0.015 in the lower layer. We assume the nearest neigh- 
bor intralayer hopping energy t = 1 and interlay hopping energy 
t± — 0.14 respectively. Only the edge-states on the right bound- 
ary are shown (red curves), (c) Band structure when the on-site en- 
ergies of the outermost column Rd in the lower layer are tuned to 
U{Rd) ~ 1. (d) Band structure when the on-site energies of the 
outermost columns in both layers are U (Rd) = U{Ru) = 1. 



a bias A between the two layers, a bulk gap can be induced 
and it has been shown that, in the insulating states, the two 
valleys carry the topolog ical charge of T2Sgn(A) in the high- 
est valance band |Ta, 17, 34]. By the above topological argu- 
ment, similar edge-states behaviors to those in the single layer 
shall be expected. In Fig.|4l we show the edge-states spectrum 
when the on-site energy is universally —0.015 in the top layer 
and 0.015 in the bottom layer. As expected, each boundary 
now carry two bands |23]. Edge band 1 has a flat dispersion 
in the region ky E [|7r, |7r], similar to the edge-band in the 



single layer case, and we find that the states near ky — n are 
completely localized on the outermost column (i?«) of the up- 
per layer Edge band 2 also has a flat dispersion near ky = tt, 
but acquires finite valley-dependent velocity near the Dirac 
points |23]. The two edge-bands have an energy difference 
equal to A = —0.03 in the flat region, which is due to the 
fact that, in edge-band 2, the states near ky = tt are localized 
on the outermost column (Rd) in the lower layer Thereby, 
the dispersion of the two edge-bands can be individually con- 
trolled by tuning on-site energies of Ru and Rd columns re- 
spectively. In particular, depending on the boundary potential, 
we can have either one or two valley-dependent gapless chiral 
modes per valley per boundary, as shown in Fig.Uc) and (d). 
In fact, the gapless edge-states when U{Ru) — U{Rd) — 1 
are equivalent to the symmetric modes confined by a line dis- 
location so that the bilayer is of AB stacking to the left of 
this topological Une defect and is of BA stacking to the right 
of the defect fsJl- The topological charge of these gapless 
modes v = —2tz in each valley, which is found consistent 
with the bulk topological charge — — r^. Indeed, in bi- 
layer graphene, the behaviors of the edge-states are also de- 
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